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Abstract. Genetic marker studies can assist restoration practice through selection of seed sources that conserve 
historical levels of gene diversity and population genetic differentiation. We examined genetic variation and structure 
within and among mainland and island populations of Eiymus glaucus, a perennial bunchgrass species native to west- 
ern North American grasslands that is targeted for grassland restoration. Island populations of E. glaucus represent 
sensitive sites and potentially distinctive seed sources for reintroduction, and little is known of their genetic compos- 
ition. Genetic diversity and structure were estimated using amplified fragment length polymorphism markers for 21 
populations and 416 individuals distributed across two coastal California mainland locations and three California 
Channel Islands. Eight primer combinations resulted in 166 markers, of which 165 (99.4 %) were polymorphic. The 
number of polymorphic bands was significantly greater among mainland populations relative to island sites, and 
locally common alleles were present for each sampled island and mainland location. Population structure was high 
(62.9 %), with most variation (55.8 %) distributed among populations, 7.1 % between mainland and island locations, 
and the remainder (37.1 %) within populations. Isolation by distance was only apparent among islands. Using marker 
data to recommend appropriate seed sources for restoration, E. glaucus seeds are best derived within islands with 
collections representing a large number of individuals from matching environments. Given the limited gene flow 
and prior evidence of adaptive divergence among populations of this species, regional collections are recommended 
in all cases to maintain diversity and to avoid long-distance introductions of highly differentiated plant material. 
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Introduction 



rhizomatous perennial exotic grasses (Mack 1989; Hayes 
and Holl 2003). California grasslands represent an ecosys- 
tem where plant community conversion is nearly com- 
plete (Mensing 1998). Mediterranean annual grasses 
and forbs introduced over the past three centuries now 
dominate the landscape, and native species exist as rem- 
nant populations in a matrix of exotics. Efforts are 
ongoing to control invasive species and to re-establish 



Widespread anthropogenic disturbance and introduc- 
tions of invasive species have resulted in the fragmenta- 
tion and conversion of grassland ecosystems worldwide 
(D'Antonio and Vitousek 1992). Temperate grasslands ori- 
ginally dominated by caespitose (bunchgrass) species 
have been altered by introduced livestock, and are highly 
vulnerable to invasion by competitive annual and 
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native perennial bunchgrasses (e.g. Moyes et a/. 2005; 
DiTomaso et al. 2007; Cox and Allen 2011), and guidelines 
are needed for the restoration of sustainable and diverse 
plant populations. 

Fragmented populations of plant species are sus- 
ceptible to environmental, demographic and genetic 
stochasticity (Shaffer 1981; Lande 1988). Restoration pro- 
grammes commonly mitigate environmental and demo- 
graphic concerns for all target species by increasing the 
number of both individuals and populations to minimize 
the probability of extirpation. In contrast, attempts to 
mitigate the loss of genetic variation and to minimize 
inbreeding have largely focused on rare and endangered 
species (Lande 1988), although interest in the genetic 
consequences of establishing or restoring populations 
of common species has grown over the last decade (e.g. 
Hufford and Mazer 2003; McKay et al. 2005; Bischoffet al. 
2010). Evidence suggests that common species are sub- 
ject to genetic erosion resulting from habitat fragmenta- 
tion at similar or even greater rates than rare species 
(Honnay and Jacquemyn 2007). Genetic variation is the 
basis of adaptation, and the loss of genetic diversity, par- 
ticularly for fitness-related traits, will impact population 
persistence as well as limit the ability of a population to 
adapt to changing environments (Frankham et al. 2002; 
Reed and Frankham 2003). Knowledge of species-level 
patterns of genetic diversity can inform and improve res- 
toration protocols when population reintroduction is a 
restoration objective. 

Primary genetic concerns for reintroduction include the 
maintenance of patterns of diversity within and among 
populations, and the introduction of genotypes adapted 
to environmental conditions at the restoration site (e.g. 
McKay et al. 2005). Population genetic structure is a func- 
tion of a species' mating system and provides an estimate 
of historical levels of gene flow and connectivity among 
locations (Slatkin 1987). Data that provide measures of 
the partitioning of genetic variation, however, do not ex- 
plain underlying causes of divergence, which may be a 
function of selection or random processes (Heywood 
1991). Direct evaluations of genotypic adaptation and 
traits under selection can be determined in common gar- 
den and reciprocal transplant studies (Linhart and Grant 
1996; Kawecki and Ebert 2004). When available, these 
studies provide valuable information for the identification 
of appropriate seed sources for reintroduction. However, 
data for common gardens are not available for all species, 
are limited in scale and may not detect all components of 
local adaptation (Nuismer and Gandon 2008). Thus, base- 
line information to describe patterns of genetic variation 
within and among populations remains relevant for res- 
toration and conservation planning, and can serve as a 



first step for management of genetic diversity in species 
reintroduction or augmentation programmes. 

The genetic consequences of seed introduction during 
restoration may have greater impacts for populations 
occupying small or isolated islands relative to mainland 
sites. Island populations often harbour lower levels of 
gene diversity and higher levels of differentiation when 
compared with the mainland, and are at increased risk 
of extinction — possibly due to greater environmental 
and demographic stochasticity (Frankham 1997). Islands 
are also disproportionately vulnerable to biological inva- 
sion, and introduced species are reported to outnumber 
native species in grasslands of the California Channel 
Islands (Halvorson 1992; Schoenherr et al. 1999; Moody 
2000). Elymus glaucus (blue wildrye) is a native bunch- 
grass species once common in California mainland 
and island grasslands (Holstein 2001; Bartolome et al. 
2004). Over the past few decades, E. glaucus has been 
a target of restoration programmes due to its wide distri- 
bution, wildlife habitat value and dense root system, 
which prevents erosion in degraded landscapes (Knapp 
and Rice 1996; Erickson et al. 2004). Two previous 
reciprocal transplant studies found evidence of ecotypic 
variation among populations of this species as a result 
of adaptation to local environments over scales of 
50-190 km (Hufford et al. 2008; Knapp and Rice 2011). 
These studies were conducted at mainland locations, 
and no information is yet available to describe genetic 
differentiation or adaptive variation for populations on 
the Channel Islands. 

In this study, we used amplified fragment length poly- 
morphisms (AFLPs) to characterize genetic diversity and 
structure among mainland and island populations of 
F. glaucus to address the following questions relevant 
for grassland restoration: (i) How much genetic variation 
is present in island populations, and how does this com- 
pare to mainland populations? (ii) Is genetic differenti- 
ation between island populations greater relative to 
nearby mainland locations? (iii) Is genetic distance corre- 
lated with geographic distance within and among islands 
and the mainland? Lower levels of diversity and strong 
genetic differentiation at island sites may indicate the 
reduced ability of island plants to adapt to altered envir- 
onments and a greater risk of local extinction of Channel 
Island populations. At the same time, the geographic 
scale of genetic differentiation serves as an indicator of 
the historical rates of gene flow. These data can assist 
with seed provenance selection for the restoration of 
E. glaucus growing in California coastal and island grass- 
lands by means of maintaining population genetic diver- 
sity and lowering the risk of introducing maladapted 
genotypes. 
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Methods 

Study species 

The genus Elymus includes 150 species distributed in 
temperate regions worldwide (Hickman 1993). Elymus 
glaucus is a perennial, non-rhizomatous bunchgrass 
with a broad geographic distribution from Canada to Mex- 
ico, and can be found throughout the western United 
States (Hickman 1993). Herbarium records for this species 
include much of the state of California (http://www. 
calflora.org), but extant populations are highly fragmen- 
ted as a result of widespread land development and 
biological invasion (Barry et al. 2006). Populations of 
F. glaucus occur in diverse habitats and plant communi- 
ties, and exhibit morphological variation across their 
range (Snyder 1951; Wilson et al. 2001). Polyploidy is 
common in the Poaceae, and E. glaucus is an allotetra- 
ploid (2n = 28) derived from Hordeum and Pseudoroeg- 
neria ancestors (Jensen et al. 1990). Previous studies of 
E. glaucus suggest that it is frequently self-pollinated 
and has a mixed mating system, allowing for some out- 
crossed pollination (Knapp and Rice 1996; Ie 2000; Wilson 
et al. 2001). Inflorescences are distinctive, narrow spikes, 
and seed dispersal is typically passive. 

Study sites and collections 

The Southern Channel Islands represent an archipelago 
of eight continental islands located at distances ranging 
from 13 to 61 km off the coast of mainland California 
(Schoenherr et al. 1999). Despite their proximity to the 
continental landmass, high numbers of endemic species 
(up to 47 % of native vegetation) are found in the island 
chain (Junak et al. 1995; Moody 2000). Introduced plants 
and animals threaten native species, and efforts to con- 
serve and restore island ecosystems are ongoing (Halvor- 
son 1994). In the present study, sampling sites included 
21 populations distributed among two mainland loca- 
tions and the three Channel Islands where E. glaucus is 
known to occur (Junak et al. 1997). The mainland sites 
were located at the University of California Sedgwick 
Reserve (Sedgwick) in Santa Ynez, California and Vanden- 
berg Air Force Base (VAFB) in Lompoc, California. Offshore 
sites were located on Santa Rosa and Santa Cruz Islands 
in Channel Islands National Park, and Santa Catalina 
Island (Fig. 1). 

Sampling locations represented diverse habitats and 
serpentine rock outcrops are common at the two main- 
land locations and Santa Catalina Island. Populations of 
E. glaucus sampled in this study occurred in oak woodland 
savannahs at Sedgwick and in open, coastal grasslands at 
VAFB. Island populations were distributed among coastal 
prairies and ephemeral, riparian environments (Fig. 2). The 
regional climate is Mediterranean, but due to the marine 



influence, coastal areas of VAFB and the Channel Islands 
experience cooler temperatures and higher humidity 
when compared with the interior of Santa Cruz Island 
and Sedgwick Ranch. Sites located at Sedgwick and in 
Santa Cruz Island's Central Valley are subject to greater 
temperature fluctuations than coastal areas, and often 
record seasonal temperature differences of 5 °C or more 
relative to the shoreline (Schoenherr et al. 1999). 

Populations of £. glaucus were identified and georefer- 
enced during several trips to each location in the spring 
and summer of 2002 and 2003 (Table 1). At each location, 
we sampled populations that represented geographically 
distant sites although some areas were inaccessible, lim- 
iting our sampling range. Efforts were made to sample 
plants separated by 0.5 m or greater within populations 
along 10- to 20-m walking transects. The patchy nature 
and small size of many E. glaucus populations, however, 
restricted the area within which we were able to collect 
leaf material, and in some cases we modified transect 
sampling to collect in a smaller radius while maintaining 
the separation of sampled plants. In general, population 
sizes ranged between 20 and 60 visible plants. At each 
site, leaves from 20 individuals were collected and stored 
in sealed bags containing silica gel for preservation. One 
population located on Santa Cruz Island was very small, 
and only allowed for 16 individual collections. 

AFLP genotyping 

Leaves stored in silica gel were transported to the Univer- 
sity of California in Santa Barbara and stored at an aver- 
age room temperature of 20 °C prior to AFLP genotyping. 
For each sampled plant, total genomic DNA was extracted 
from ~20 mg of silica-dried leaf tissue using the DNeasy 
Plant Mini Kit (Qiagen). A total of 416 plants representing 
21 £. glaucus populations distributed among 9 mainland 
and 12 island sites were included in subsequent marker 
analyses. 

Molecular markers were amplified following the proto- 
col of Vos et al. (1995) with little modification. Approxi- 
mately 250 ng of DNA were digested with EcoRI and 
Msel restriction enzymes and ligated to corresponding 
adapters. The restriction -ligation mix was diluted 1 : 10 
and polymerase chain reaction (PCR)-amplified with 
EcoRI-A and Msel-C pre-selective primers. The resulting 
PCR template was diluted 1:10 prior to amplification 
with selective EcoRI primers that were 5' end labelled 
with IRDye 700 or 800. Eight primer combinations were 
selected due to clarity and repeatability of bands 
(Table 2). Amplification products of duplexed selective 
PCR reactions were denatured and separated on 7 % 
acrylamide gels using a LI-COR 4200 DNA sequencer 
(LI-COR, Lincoln, NE, USA). One or more duplicate samples 
were routinely included in gel runs for quality control. The 
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presence or absence of AFLP bands was scored manually 
using SAGA 2.0 MX software, and ambiguous bands were 
recorded as missing. Error rates (calculated as the num- 
ber of mismatched genotypes divided by the number of 
bands) were 3-4 %. 

Genetic analysis 

We analysed AFLP marker data using methods previously 
employed for allotetraploid species, including band- 
based and allele-frequency metrics (e.g. Kang et al. 
2005; Wagner et al. 2012). Band-based metrics compute 
distance-based measures of similarity within and among 
populations from the matrix of marker presence/absence 
data without inferring allele frequencies (Bonin et al. 
2007). In contrast, allele-frequency methods estimate 
standard population genetic statistics; these methods 
were developed for diploid organisms and may not be 
valid for polyploid species. Allele-frequency methods 
were tested by Wagner et al. (2012) for an allotetraploid 
species, and estimated values correlated strongly with 



band-based metrics, since allotetraploid species likely 
undergo meiotic segregation within each parental gen- 
ome similar to diploid species (Soltis and Soltis 1993). 

Genetic diversity 

Using a band-based approach, we calculated the number 
and proportion of polymorphic bands (PPBs) within and 
among island and mainland sites by means of FAMD 
1.2 software (Schluter and Harris 2006). We also com- 
puted the number of fixed (or monomorphic) bands. 
The number of private and locally common bands 
(restricted to a limited area and found in <25-50 % of 
populations) was determined in GenAlEx 6.5 (Peakall 
and Smouse 2006). Differences in the number of poly- 
morphic or fixed bands observed for island and mainland 
sites were compared using generalized linear models 
and a Poisson distribution in JMP 9 statistical software 
(SAS Institute Inc., Cary, NC, USA). The proportion of 
polymorphic loci (PLP) and expected heterozygosity (H e ) 
were estimated for comparison to band-based metrics 
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Figure 2. Riparian site representing potential E. glaucus habitat on 
Santa Rosa Island. 

using a fragment-frequency approach in AFLP-SURV 1.0 
(Vekemans 2002). This method assumes fixed homozy- 
gosity resulting from self-pollination and may overesti- 
mate marker frequencies for outcrossing taxa, but 
would meet expectations for F. glaucus. We subsequently 
compared band-based metrics with allele-frequency es- 
timates using regression analysis. All files for genetic 
analyses were prepared using functions available in 
GenAlEx or the package AFLPdat (Ehrich 2006) for R soft- 
ware (R Development Core Team 2011). 

Genetic differentiation among populations 
and regions 

We conducted a hierarchical analysis of molecular vari- 
ance (AMOVA) in GenAlEx to partition genetic variation 
within and among populations nested in two regions: 
island and mainland. AMOVA components of variance in- 
clude <J> PT , which is considered to be an analogue of F S t 
(Wright 1951; Peakall and Smouse 2006). The significance 
of the different components of variance was tested with 
9999 permutations. Values of 4> PT were also calculated 
separately for populations within island or mainland re- 
gions. A Mantel test (10 000 permutations) was used to 
determine if there was an association between genetic 



distance measured as a matrix of linearized <S PT values, 
and logio transformed geographic distances (Rousset 
1997). The Mantel test was also conducted for the 
two subsets of island and mainland data. The indirect 
rate of gene flow (N m ) was estimated following Wright 
(1951), and pairwise genetic distance (4> PT ) values were 
calculated and assayed for significance with 9999 
permutations. 

We compared AMOVA results for population genetic 
structure with allele-frequency estimates calculated in 
HICKORY software 1.1 (Holsinger et al. 2002; Holsinger 
and Lewis 2003). HICKORY employs a Bayesian approach 
to estimate (comparable to F S Tand 4>p T ) using domin- 
ant markers, and does not assume Hardy- Weinberg 
equilibrium. We computed for three alternative mod- 
els (full model, f = 0 model and f-free model) and 
selected a suitable model with the lowest deviance infor- 
mation criterion (DIC). Markov chain Monte Carlo (MCMC) 
parameters were set to default values (burn-in 50 000, 
sampling 250 000). To test for significant differences in 
population genetic structure among island and mainland 
locations, we ran posterior comparisons of 6 iu) values. 
If 95 % confidence intervals for the difference in paired 
samples included zero, we concluded that values 
were not significantly different. 

To test the assumption that the 21 F. glaucus sites re- 
presented distinct populations, we used Bayesian cluster- 
ing methods implemented in STRUCTURE 2.3.3 software 
to assign individuals to populations by employing the re- 
cessive alleles option for dominant markers (Pritchard 
et al. 2000; Falush et al. 2007). We ran 20 iterations of 
each K=l-23 possible clusters using the default 
model that infers alpha and assumes admixture and cor- 
related allele frequencies. Every run included a burn-in 
period of 150 000 MCMC cycles and 300 000 MCMC itera- 
tions (University of Oslo Bioportal; Kumar etal. 2009). The 
most likely number of clusters represented by the AFLP 
data was identified using the method described in 
Evanno et al. (2005) and implemented in HARVESTER soft- 
ware (Earl and vonHoldt 2012), which calculates AK as the 
second-order rate of change of the log probability of the 
data. CLUMPP 1.1.2 software (Jakobsson and Rosenberg 
2007) aligned the 20 replicate runs and results were plot- 
ted with DISTRUCT 1.1 software (Rosenberg 2004). 

Results 

Genetic diversity 

We scored clear and unambiguous AFLP bands as present 
or absent for each sampled individual. Of the 166 AFLP 
markers, 165 (99.4 %) were polymorphic for the full data- 
set. The proportion of missing data was calculated at 
1.91 %. The average number of scorable bands generated 
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Table 1. Sampled locations and genetic diversity indices for 21 E. glaucus populations including latitude (N ) and longitude (W ), sample size (n), 
number of locally common (f ) and fixed bands (FB), per cent polymorphic bands (PPB), per cent polymorphic loci (PLP) and expected 
heterozygosity (H e ) with standard errors 



Location 


N° 


W° 


Description 


Population ID 


ft 


f 
j 


FB 


PPB (%) 


PLP (%) 


Mean 


SE 


jUIILU luluim iu 


33 20.81 


1 1 O lD.Dj 


RiillnicH f nni/nn 
DUUIUbll L-Uliyuil 


CI 


20 


16 


75 


32.7 


30.9 


0.120 


0 015 




33 21.43 


118 25 57 


LU|JC L.UI iyui 1 


02 


20 


15 


59 


50.9 


50.6 


0 118 


0.012 




33 20.85 


118 24.32 


MiHHIo f nn\/nn 
l v ll(JU.lc l_UIIVUII 


rj3 


20 




31 


72.7 


72.3 


0.213 


0 015 




33 21.19 


118 21 63 


nuy ui ebb 


C4 


20 


13 


40 


60.6 


56.0 


0 157 


0.013 


JUI 1 LU 1 U L 


34 00 75 


119 47 77 


DnrtoTi lain 
rUI LciUclU 


SCI 


20 


8 


82 


25.5 


25.3 


0.095 


0 014 




34 00.35 


119 44 96 


\/nllo\/ Pnnrl 
VUllcV l\UUU 


SC2 


20 


4 


97 


18.8 


18 7 


0.075 


0.013 




34 00 94 


119 35 55 


Fnn nf mo 1 Ino 
CIIU UI Lilt; LMIC 


SC3 


20 


2 


77 


28.5 


28.3 


0 101 


0 014 




34 02 44 


119 34 45 


JL.UIUIUII V_UIIVUII 


SC4 


16 


8 


100 




4.8 


0.013 


0.005 


jUIILU r\UbU 


33 59.82 


120 05.25 


1 nhn f nmmn 
LUUU i_uiiyuii 


SRI 


20 


4 


80 


23.6 


22.4 


0.053 


0.009 




3 ~i Era Q"7 
55 Dy.o/ 


1 ir\ n3 ~ir\ 


Cherry Canyon 


5KL 


ZU 


5 


"7 "7 




5\i.l 


u.uoy 


U.U1U 




33 58.97 


120 02.98 


Water Canyon 


SR3 


20 


5 


58 


47.3 


47.0 


0.117 


0.013 




33 58.55 


120 01.00 


Box Canyon 


SR4 


20 


7 


38 


61.2 


59.8 


0.150 


0.013 


Sedgwick 


34 43.32 


120 02.16 


Figueroa 4 


SI 


20 


4 


79 


23.0 


22.4 


0.060 


0.011 




34 43.24 


120 02.18 


Figueroa 3 


S2 


20 


8 


55 


49.1 


48.8 


0.089 


0.009 




34 43.00 


120 02.34 


Figueroa 2 


S3 


20 


6 


52 


51.5 


51.2 


0.141 


0.013 




34 42.62 


120 02.40 


Figueroa 1 


S4 


20 


4 


97 


18.2 


18.1 


0.035 


0.008 




34 41.43 


120 02.75 


Ranch House 


S5 


20 


13 


36 


67.3 


66.9 


0.211 


0.015 


VAFB 


34 48.08 


120 30.86 


Campground 


VI 


20 


8 


61 


47.3 


47.0 


0.143 


0.015 




34 36.61 


120 31.90 


Pasture 


V2 


20 


10 


62 


50.3 


50.0 


0.123 


0.012 




34 35.64 


120 31.91 


San Miguelito 


V3 


20 


15 


31 


65.5 


65.1 


0.141 


0.012 




34 34.30 


120 30.92 


Sudden/Ave I 


V4 


20 


13 


87 


24.8 


24.7 


0.070 


0.011 










Island 


236 




67.8 


98.8 


37.2 


0.107 


0.015 










Mainland 


180 




62.2 


96.4 


43.8 


0.113 


0.018 










Mean 








40.7 


40.0 


0.109 


0.011 



Table 2. Combinations of EcoRI and Msel selective primers, and the 
total number of bands scored and per cent polymorphism (PLP + SE) 
generated by each primer combination for E glaucus populations 



Primer pairs 


# Bands 


% PLP 


% SE (PLP) 


E-AAC/M-CAG 


17 


38.38 


4.93 


E-ACC/M-CAG 


17 


37.54 


4.40 


E-AAG/M-CGG 


22 


39.39 


4.65 


E-AGC/M-CGG 


23 


43.06 


5.00 


E-ACA/M-CAG 


25 


41.33 


4.82 


E-AGG/M-CAG 


19 


37.59 


5.17 


E-AAC/M-CGG 


23 


35.40 


4.32 


E-ACA/M-CGG 


20 


46.19 


5.72 


Mean 


20.75 


39.96 


4.11 



by each AFLP primer combination was 20.75 (range of 
17-25; Table 2). The number of polymorphic bands 
declined significantly for subsets of island or mainland 
data, reducing markers for analysis in some cases by 
more than half (Table 1). Mean genetic diversity for all 
samples was relatively high (PPB = 40.7 %), but varied 
considerably among populations (6-73 %). The lowest 
reported values for polymorphism were recorded for the 
population (SC4) on Santa Cruz Island with only 16 indivi- 
duals, and may represent a recent founder event. 

All values describing genetic diversity (PPB, PLP and H e ) 
were strongly correlated (pairwise comparisons, r > 0.91, 
P< 0.0001), indicating concordance among band-based 
and allele-frequency metrics. The average expected het- 
erozygosity (H e ) within populations was low (0.1093) and 
varied among populations and locations (Table 1). The 
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average expected heterozygosity within locations was 
greatest for the Catalina, VAFB and Sedgwick populations 
(0.152, 0.119 and 0.107, respectively) and declined for 
Santa Cruz (0.071) and Santa Rosa (0.097) Islands but 
was not significantly different among mainland and 
island locations. Comparisons of genetic diversity para- 
meters between the two groups of mainland and island 
populations detected a significantly greater number of 
polymorphic bands among mainland populations (mean 
of 73 vs. 63 bands, respectively; P = 0.007). The mean 
number of fixed bands among populations did not differ 
significantly between mainland and island regions at the 
P = 0.05 level (67.8 vs. 62.2 bands; P = 0.101). Two pri- 
vate bands separated pooled island and mainland sites, 
and locally common bands were present among popula- 
tions at each of the three island and two mainland loca- 
tions (Table 1). 

Genetic differentiation among populations 
and regions 

Most of the variation in the AFLP profiles reported 
here represented variation among populations within 
regions. Hierarchical AMOVA for the island and main- 
land data partitioned 37.1 % of the variation within 
populations, 55.8 % among populations within regions 
and 7.1 % between regions (Table 3). All values were 
significantly different from zero (P< 0.0001) and 4> PT 
summed to 0.629. The average pairwise 4> PT for all 
sites was 0.614 with a range from 0.159 (VAFB popula- 
tions, V3 and V4) to 0.91 (Sedgwick and Santa Cruz 
populations). All pairwise comparisons for $ PT were 
significant [see Supporting Information]. Population 
genetic structure declined within regions (<X> PT = 0.422 
among mainland and 0.400 among island populations) 
relative to structure calculated for all sites. The results 
of the Mantel test of linearized <t> PT values, and logio 
transformed geographic distances indicated no support 
for isolation by distance among all island and mainland 
sites (P = 0.175). This result was repeated when 



Table 3. Analysis of molecular variance (AMOVA) results for E. glaucus 
populations located within and among the two (island and mainland) 
regions 

Source of d.f. Sum of Variance Variation 

variation squares components (%) 

Among regions 1 762.55 2.08 7.1 

Among populations, 19 336.02 16.42 55.8 
within regions 

Within populations 395 10.95 10.95 37.1 



mainland sites were considered alone (P = 0.293). How- 
ever, isolation by distance was detected among popula- 
tions sampled on the three Channel Islands (matrix 
correlation coefficient = 0.153, P = 0.023). This result 
was limited to the groups of populations among islands 
and isolation by distance was not observed for popula- 
tions within each island [see Supporting Information]. 
The value of N m , or the indirect rate of gene flow, was 
calculated as 0.147, indicating that the number of mi- 
grants each generation was significantly <1 and that 
genetic drift plays a role in population differentiation 
(Slatkin 1987). 

Using the Bayesian approach, two models of 0 lu) (full 
and f = 0) resulted in low DIC values, and we selected 
the full model with the expectation that the inbreeding 
coefficient (f) is >0 in this autogamous species (K. Hol- 
singer, pers. comm.; Wilson et al. 2001). Under the full 
model, the value of for all 21 sites was 0.563, 
while values for the subsets of mainland or island sites 
were 0.557 and 0.556, respectively. Comparisons of pos- 
terior distributions detected no significant difference in 
population structure between island and mainland 
regions (difference in values: 0.001; confidence 
interval: -0.0379, 0.03932). Estimates of 0 (II> were not 
greater for populations located on islands relative to 
sampled populations on the mainland. Sedgwick and 
VAFB also did not differ significantly for estimates of 
However, Santa Catalina and Santa Rosa popula- 
tions were significantly less differentiated than popula- 
tions from Santa Cruz (difference in 9 iu) values: 0.189; 
confidence interval: 1349, 0.2445). Genetic structure 
was greater among sampled sites located on Santa 
Cruz Island. 

Not all sampled populations were genetically distinct 
from one another. The highest likelihood of the number 
of genetic clusters represented in the AFLP dataset was 
consistently obtained with K = 16 (Fig. 3). In effect, the 
21 collections represented 16 genetically distinct popu- 
lations. Sites V3 and V4 corresponded to a single genetic 
cluster and these sites were located ~3 km apart. Indi- 
viduals among the four collection sites at Santa Rosa 
Island were also highly similar (2-7 km apart) despite 
their location in separate drainages on the northeastern 
shore of the island. Lastly, pairs of sites within Sedgwick 
and within Santa Cruz Island were overlapping. 
Altogether, evidence of admixture was low and plants 
derived from no more than five to seven of the original 
sites appeared to have patterns of diversity representa- 
tive of dispersal and gene flow within islands (Fig. 3). 
Essentially, these individuals shared membership in 
more than one population. Most populations had 
fixed genetic differences that did not vary widely 
among individuals. 
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Discussion 

Strong levels of genetic differentiation observed in this 
study indicate that £. glaucus populations are highly 
self-pollinating and isolated at small distances. As a re- 
sult, populations represent genetically homozygous 
lines and the distribution of diversity is heavily influ- 
enced by the breeding system (Schemske and Lande 
1985). Similar patterns have been observed for previous 
studies of E. glaucus in North America, as well as among 
populations of other grass species worldwide (e.g. Jain 
and Allard 1960; Nevo et al. 1986; Knapp and Rice 
1996). Predominantly self-pollinating species lead to 
different expectations for patterns of genetic diversity, 
and represent unique challenges for conservation that 
merit consideration when these species are targeted 
for reintroduction. In our study of E. glaucus among Cali- 
fornia Channel Island and mainland locations, we ad- 
dressed questions regarding the levels and patterns of 
genetic diversity and differentiation at relatively iso- 
lated island sites, and their consequences for ecological 
restoration. 



Genetic diversity 

Despite strong evidence for self-pollination, marker data 
indicated that considerable genetic variation is present in 
blue wildrye. Measures of the mean proportion of poly- 
morphism within studied £. glaucus populations were 
~40 %. The distribution of genetic diversity may be the 
result of factors other than inbreeding, including the poly- 
ploid origin of £. glaucus as well as the potential, however 
small, for gene exchange among populations (Jain and 
Allard 1960). Genetic polymorphism present in alloploid 
grasses such as £. glaucus may be a consequence of seg- 
regation between homologous chromosomes within 
each parental genome, resulting in fixed heterozygosity 
in selfing lineages (Soltis and Soltis 1993). Additionally, 
self-pollination is rarely complete and gene exchange 
among lineages is possible through wind pollination in 
combination with seed dispersal among populations. Pre- 
vious isozyme studies of £. glaucus growing in the north- 
western United States and British Columbia also detected 
high levels of polymorphism at the species (77-80%) 
and population level (22-31 %) (Knapp and Rice 1996; 
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Ie 2000; Wilson et al. 2001), suggesting that high 
levels of genetic variation are maintained across this 
species' range. 

Relative to the entire region sampled in this study, 
measures of genetic variation declined significantly 
when data were partitioned within island or mainland lo- 
cations and particularly for individual populations, many 
of which were genetically homogeneous for AFLP mar- 
kers. Mean expected heterozygosity (H e = 0.109) for 
Channel Island and coastal California populations was 
low and corresponded to values recorded for selfing 
species with passively dispersed seeds (H e = 0.097; 
Hamrick and Godt 1996), but varied among sampled 
sites. Populations located on Santa Cruz and Santa Rosa 
Islands had the lowest average levels of polymorphism 
while populations on Catalina Island were the most di- 
verse. Levels of genetic diversity may vary among island 
and mainland locations, perhaps as a result of variation 
in the breeding system in addition to stochastic demo- 
graphic and environmental factors (Schemske and 
Lande 1985). 

Comparisons of the number of polymorphic bands be- 
tween the two regions representing islands and the main- 
land indicated that genetic diversity was significantly 
higher for mainland populations relative to the Channel 
Islands. These data support predictions that island popu- 
lations are genetically depauperate as a result of founder 
events, low rates of gene flow and subsequent genetic 
drift (Frankham 1997). Sampled populations for this 
study commonly represented fewer than 60 visible plants 
and locating populations on the islands was difficult, 
probably as a result of greater geographic isolation 
among £. glaucus populations on the islands. The conse- 
quences of relatively lower levels of diversity remain un- 
clear, however. In many taxa, genetic diversity strongly 
correlates with population fitness and may be useful for 
identification of populations at risk of decline (Reed and 
Frankham 2003), but this pattern requires further investi- 
gation in self-pollinating species (Leimu et al. 2006). High 
levels of inbreeding can purge deleterious alleles so that 
selfing species are unlikely to experience inbreeding de- 
pression, and yet heterosis is commonly reported when 
gene exchange occurs (Jain and Allard 1960; Schemske 
and Lande 1985). 

Genetic differentiation among populations 
and regions 

We discovered strong genetic differentiation in E. glaucus 
within and among island and mainland sites. Most of 
the variation was distributed among populations within 
regions and overall estimates of F S j were very high 
(d> PT = 0.629, = 0.563) and similar to all prior studies 
of this species (Knapp and Rice 1996; Ie 2000; Wilson 



et al. 2001). These data point to the dominant role of 
the breeding system in shaping the population structure 
of £. glaucus. At the same time, relative levels of isolation 
may also affect patterns of genetic differentiation if the 
distance among populations further limits rare episodes 
of gene flow. 

Tests for isolation-by-distance (IBD) were not signifi- 
cant for populations on the mainland and IBD was only 
apparent between islands. Two previous mainland stud- 
ies of £. glaucus conducted at a larger scale also failed 
to detect IBD (Knapp and Rice 1996; Wilson et al. 2001), 
while low levels of IBD were detected for sites sampled 
throughout British Columbia (Ie 2000). Overall, genetic 
structure was not greater among Channel Island popula- 
tions than among mainland populations and IBD was not 
observed among populations within each island. 

Bayesian cluster analysis presented some evidence for 
admixture in E. glaucus as the result of gene flow, al- 
though this evidence was limited to few individuals and 
populations (e.g. SR4, Santa Rosa Island). Inference of 
genetic structure detected 16 genetically distinct popula- 
tions among 21 collections, and confirmed strong popu- 
lation differentiation in £. glaucus among islands and the 
mainland. The average pairwise geographic distance 
among sites was ~100 km. However, many sites within 
locations were only 2-3 km apart. We noted that the 
two closest collection sites (SI and S2; 0.15 km apart) 
were strongly genetically differentiated. High levels of dif- 
ferentiation at small distances are common among self- 
pollinating grass species and represent challenges for 
their restoration in degraded ecosystems (Allard 1975). 

Implications for restoration 

Evidence for lower levels of polymorphism among islands 
suggests that island populations of blue wildrye are gen- 
etically depauperate relative to the mainland. In these 
cases, the introduction of diverse seed sources in restor- 
ation has been proposed to increase population viability 
by reduction of inbreeding depression, and creation of 
opportunities for hybrid vigour (Broadhurst et al. 2008; 
Pekkala et al. 2012). Conversely, local seed collections 
are often proposed to avoid disruption of local adaptation 
and the incidence of outbreeding depression among hy- 
brid progeny of remnant populations and introduced 
plants (Hufford and Mazer 2003). How do these contrast- 
ing views apply to restoration of £. glaucus? 

If inbreeding depression is weak in self-fertilizing 
species (Schemske and Lande 1985), the preservation of 
genetic variation to conserve population viability should 
no longer be a large concern in reintroduction pro- 
grammes. It does not follow, however, that the diversity 
of lines represented in seed mixes would also be of little 
concern. There is strong evidence for the maintenance of 
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heterozygotes in selfmg species, and heterozygote ad- 
vantage has been proposed as a mechanism to improve 
population fitness when outcrossing is rare (Allard and 
Jain 1962). Ideally, when reintroducing species such as 
blue wildrye, the goal would be to maintain or even in- 
crease genetic variability among plants. In this case, 
the introduction of seeds representing large numbers of 
mixed lines would be supported (Broadhurst et al. 
2008). However, non-local genotypes may represent 
plants poorly adapted to local environmental conditions. 

High levels of genetic differentiation may be the prod- 
uct of both genetic drift and natural selection (Slatkin 
1987), and evidence for local adaptation in E. glaucus 
has been detected in two reciprocal transplant studies 
(Hufford et al. 2008; Knapp and Rice 2011). In addition, 
other studies of blue wildrye detected high levels of 
fixed variation that were correlated with morphological 
variation among sites (e.g. Knapp and Rice 1996; Wilson 
et al. 2001; Erickson et al. 2004). In light of this evidence, 
population genetic divergence apparent in E. glaucus may 
reflect not only limited gene flow but also adaptation to 
heterogeneous environments. We noted that genetic 
structure declined when the full dataset was partitioned 
among islands and mainland locations (4>p T = 0.422 for 
the mainland and 0.400 for island populations relative 
to 0.629 among all sampled populations), and this scale 
corresponds to the scale of adaptive differentiation de- 
tected in field studies of this species. We consequently 
recommend the use of multiple seed collections specific 
to each Channel Island (for use at that island) to maintain 
diversity while avoiding long-distance introductions of 
non-local genotypes. Given low levels of pollen flow, 
these introductions are not likely to affect the viability 
of existing populations, and in rare cases where gene 
flow does occur, hybrid vigour among progeny would 
likely result in the formation of new selfmg lines (Allard 
and Jain 1962; Nevo et al. 1986). Similar guidelines 
might apply to mainland locations, but would improve 
with further investigation of mainland sites. 

We finish by noting that our recommendations are 
based on both evidence and speculation. It is unlikely 
that AFLP markers, which are presumed neutral, will cor- 
respond to adaptive differentiation in this species, and 
strategies for transfer of E. glaucus are derived from pre- 
vious studies of adaptive variation as well as marker data 
(Kirk and Freeland 2011). Delineation of seed transfer 
zones, or regions within which translocation of plant 
materials is unlikely to result in the introduction of mal- 
adapted genotypes, is currently the best method to pre- 
dict adaptive divergence and select suitable germplasm 
for restoration (e.g. Johnson et al. 2004). Yet, data for 
adaptive differentiation at the scale of contemporary 
grassland restoration are scarce. Additional research 



(e.g. Hufford et al. 2008) to examine the scale of adaptive 
divergence among populations in sensitive regions would 
benefit restoration of temperate grasslands. In the interim, 
marker data can assist with the selection of germplasm 
to minimize the risk of introductions and yet maintain 
high levels of genetic diversity. Restoration sites, includ- 
ing large-scale seeding efforts under way across the 
western United States (Burton and Burton 2002), will 
serve as the proving ground for marker predictions of 
short- and long-term population viability in contempor- 
ary and changing environments. 
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